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ABSTRACT 


The cross sections (CS) for capture collisions between the symmetric top molecule 
CH 3 CN and molecular ions are calculated numerically. These ion -molecule collisions 
have been studied experimentally and have large reaction CS, Q R . Calculations are done 
to study the role of the added degree of rotational freedom by comparing capture CS for 
tops with linear molecules. The monoenergetic CS CT c (^j) (where €j is the ion energy 
and T r is the top rotational temperature) are less than the maximum classical 
CS ^ max > but much larger than Langevin CS The integrated a , Q c< agree well with 
Q r for the collisions studied. The o Q values for CH^CN are nearly independent of 
target geometry. 
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CALCULATION OF CAPTURE CROSS SECTIONS FOR ION- POLAR- 
MOLECULE COLLISIONS INVOLVING METHYL CYANIDE* 
by John V. Dugan, Jr., James H. Rice, and John L. Magee* 
Lewis Research Center 

SUMMARY 


The cross sections for capture collisions between the symmetric top molecule 
methyl cyanide (CHgCN) and several molecular ions are calculated by a numerical 
approach previously used to study collisions between ions and linear polar molecules. 
These ion -molecule collisions involving CHgCN have been studied in the mass spectrom- 
eter and have large experimental reaction cross sections. Calculations are also done 
for symmetric tops with varying moments of inertia about their symmetry axes. This 
was done to study the role of the added degree of rotational freedom by comparing 
symmetric top capture cross sections with those for linear molecules. The mono- 
energetic cross sections obtained in this report (for tops at rotational temperature 
Tj^) are less than the maximum classical cross section a max which is assumed by the 
experimentalists. These a values are much larger than Langevin cross sections cr T 
and have a slope ej ‘ 03 where is the ion energy. The integrated values of a c 
(macroscopic cross section Q c ) are within 50 percent of the reaction cross sections for 
the collisions of CHgCN + and CgD^ 2 + with CHgCN. The cross sections for different 
symmetric tops are nearly equal to the linear molecule results so that a c appears 
independent of molecular geometry of the target species. 


INTRODUCTION 

The capture cross sections for collisions between ions and polar molecules have been 

A summary of the results of this report was presented at the Annual Meeting of 
the American Physical Society, Chicago, 111., Jan. 29-Feb. 1, 1968. 

^Professor of Chemistry, University of Notre Dame. 



studied by numerical methods in references 1 to 3. In reference 1 the species studied 
were hypothetical linear and symmetric top molecules. The accuracy of the numerical 
integration technique was checked with the conservation of energy equation. An auto- 
matic step -size changing routine was used to calculate the capture cross section for 
theHCl + + HC1 collision in reference 3. The numerical results (plotted against relative 
translational energy) showed the same slope as an approximate theory predicts, but the 
absolute values differed by about 40 percent. The ion-molecule (IM) trajectories were 
studied under time reversal to ensure correct solutions of the differential equations. 

The rotational motion of linear HC1 was followed in such a way as to avoid numerical 
singularities that arise from retaining only two coordinates in any integration step 
(ref. 3). 

Certain ion-symmetric top (polar) molecule reaction cross sections have been 
studied in the mass spectrometer as reported in references 4 and 5. These reaction 
cross sections can be related to the capture cross sections by reaction efficiencies. The 
methyl cyanide molecule (CHgCN) is a polar target species which was particularly well 
studied in reference 5. The purpose of this paper is to calculate numerical capture 
cross sections for certain IM collisions involving the CHgCN molecule. 

It also includes a description of the variable step -size approach, used both in this 
study and in reference 3, with details of its application. The cross sections are 
given as functions of ion translational energy for a distribution of targets at rota- 
tional temperature Tj^. Integrated values of ct c will be compared with experimental 
data, although the comparison is somewhat qualified. The computer approach employs 
the integration routine of reference 3 with a special Euler angle transformation. This 
coordinate transformation ensures accurate solutions for the symmetric top rotational 
motion. 


DETAILS OF CALCULATION 
Numerical Approach 

The Lagrangian equations of motion were solved for the ion -molecule (IM) collision 
with an interaction potential composed of the classical ion -permanent dipole and ion- 
induced dipole (electronic polarizability terms). All collisions studied involved CHgCN, 
which was represented as a spherically polarizable ideal permanent dipole target inter- 
acting with point-change ions. The molecule rotates as the IM pair move relative to 
each other translationally. These equations and the values of molecular constants are 
included in appendix A and table I. The coordinate system used for the collision model 
is shown in figure 1. The numerical singularity difficulties (whose avoidance requires 
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TABLE I. - MOLECULAR CONSTANTS 


[Dipole moment, 3.92 DU (1. 31xl0“ 29 
(C)(m)); electronic polarizability, 

3. 8’A 2 (3. 8xl0~ 29 m 2 ); moment of 
inertia, Ip 9. 12xl0 -46 (kg)(m 2 ) 
for all cases.] 


Case 

Reduced mass, 
m, 
kg 

Moment of inertia, 

h’ 

(kg)(m 2 ) 

I 

3.41xl(T 26 

5.47X10' 47 

II 

4.78 

5.47X10' 47 

m 

3.41 

2.19X10' 46 

IV 

3.41 

0 



Figure 1. - Coordinate system employed in computer study of interaction between ion and 
symmetric top polar molcule methyl cyanide. 
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the coordinate transformation) are caused by the sin £ term in the denominator of 
equation (A7). This term prevents accurate calculation of the azimuthal acceleration fj. 
The singularity is avoided by transforming to a second Euler angle system whose z-axis 
is orthogonal to that of the first system. The details of this transformation have been 
given in reference 6. The transformation was made whenever the rotating top entered 
"polar caps" within 30° of either the positive or negative z axis (in either coordinate 
system). 

The integration procedure and step-size changing routine are based on the relative 

error encountered in conservation of energy. These are described in appendix B. A 

-2 -3 

total error (in energy) in the range 10 to 10 was generally sufficient to assure satis- 
factory behavior of coordinates and derivatives under time reversal. When energy 
conservation was violated even with small step size (low -ion velocity, small -impact 
parameter) the order of the integration routine was changed to solve the trajectory 
satisfactorily. 


Experimental Reaction Cross Sections 

Experimental reaction cross sections for ion-polar molecule systems have been 
obtained from self-consistent fits of cross section against ion energy in reference 5. It 
is assumed in that study that the collision pairs have variable reaction efficiencies. It 
should be noted that the present computer calculation takes no account of such reaction 
efficiencies (which are called transmission coefficients in ref. 5). The capture cross 
section is simply a numerical upper limit to the reaction cross section (see eq. (9) 
of ref. 5). The classical upper limit to the capture cross section has been discussed 
previously in references 2 to 6. 

It should be noted that reactions can occur at any value of ion energies from a 
thermal lower limit to the maximum energy (set by the applied repeller voltage in the 
experiment). Thus the corresponding numerical cross section of interest is an inte- 
grated value of the microscopic cross section as a function of ion energy. These in- 
tegrated cross sections will be compared with experimental cross-section results. 
However, the manner of interpretation of results is limited because the numerical cal- 
culation takes no account of reaction efficiencies. This comparison is further compli- 
cated because the approximate form of the capture cross section adopted in reference 5 
is not valid at low ion energies as will be discussed later. 
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Initial Conditions 


The initial conditions for the dipole rotator were generated from a distribution at 
Tj^ = 500 K because experiments were done at this temperature. Three to five sets of 
36 and/or 54 trajectories were used to generate each point on the C 0 against b 2 

4 C ^ r 

curve. This was done for initial ion velocities of 5. 5X10 , 10 , and 2x10 centimeters 
per second. Values of reduced mass m = mjmg/(m^ + mg) (where m^ is the ion mass 
and m 2 is the mass of the polar molecule) correspond to the CH^CN + CHgCN* and 
CHgCN + CgDjg collision pairs. The ion energies corresponding to the velocities 
studied range from several hundredths of an electron volt (eV) to several eV. The IM 
reactions of interest have been studied experimentally from ion energies ~kT R up to 
several eV. 

Four different collision pairs were studied. The first two systems consisted of 
CHgCjj colliding with (I) CHgCN* (the "parent” ion) and (II) with CgDjg. Case in con- 
sisted of a pseudo CHgCN molecule with a moment of inertia about the symmetry axis 
Ig equal to four times that of the value for CHgCN. The limit where Ig goes to zero 
was studied in case IV. This system consisted of a linear polar target with the same 
polarizability and dipole moment as CHgCN. The reduced mass assigned to the pairs in 
cases HI and IV was the same as for the CHgCN molecule and ion (case I). 

Numerical Capture Cross Section 

Because the probability of capture is not a simple function of impact parameter, the 
numerical capture cross section a c must be defined in terms of a collision ratio C^ 
as in references 1 and 3. This collision ratio is the fraction of IM trajectories in 
which the ion and molecule approach within 2 A (0. 2 nm). The choice of 2 A as the cut- 
off distance r £ appears arbitrary, but in the present calculation it has been found that 
a is insensitive to r as long as the r value is not made too large. This limit on 
r c is <4 A (<0. 4 nm) at low ion velocity and less than 3 A (0.3 nm) at high ion velocity. 
These r c values correspond to the range of van der Waals radii and hard-core cross 
sections reported for the molecules of interest in reference 5. 

The numerical cross section as a function of the free -particle impact parameter b 
is given by 

a c s * b c = 1T C rW (1) 

o 

that is, simply ir times the area under the curve C R against b . For comparison 
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with experiment one requires the integrated numerical cross section (see eq. (5) of 
ref. 5). 


Q c 


1 


"m 



( 2 ) 


In equation (2), e m is the maximum ion energy (fixed by the repeller voltage), and 
is the minimum ion energy (taken as thermal). 

The maximum classical cross section for capture has the following form as men- 
tioned in references 1 to 5. 


_ njie _ 

max L D L 


(3) 


where e is the relative translational energy of the IM pair and a, is the familiar 

2 1/2 ^ 

Langevin cross section which is equal to n(2 cte /e) . The experimentalists have 

adopted this form of the cross section to obtain reaction efficiencies from their data 
The relative translational energy is related to the ion energy for a fixed velocity by 
e = (m/mj)e^. Thus <^ max can be rewritten as 


a max^ e l^ 



a L (e l> 


(4) 


where cr^^) = rcpe/ej and a^gj) = n(2ae 2 . It is assumed that the ion velocity 
is much greater than the neutral velocity; this introduces some slight error for energies 
near Also, it will later be pointed out that the term containing e i is not negligible 
in the cr^ contribution to the integrated cross section Q max - 


RESULTS AND DISCUSSION 
Case I: CH 3 CN + + CH3CN Collision 

Curves of collision ratio C 0 against impact parameter b for the CH^CN-parent 

K 4 5 0 5 . 

ion collision are given in figure 2 for ion velocities of 5. 5x10 , 10 , and 2x10 centi- 
meters per second. The variations in C R with specified total error and random num- 
ber sets (36 or 54 trajectories) are shown by the range -of -error bars. The v = 5. 5xl0 4 
value corresponds to an ion energy of only 0. 064 eV. The somewhat irregular curve in 
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Collision ratio, 



0 10 20 30 40 50 60 

Impact parameter squared, b 2 , A 2 (or 10" 2 ^ m 2 ) 

(c) Ion velocity, 2xl0 5 centimeters per second. 

Figure 2. - Variation of collision ratio with free-particle impact parameter for 
rotational temperature of 500 K. Methyl cyanide parent-ion collision; dipole 
moment, 3.92 Debye Units (1. 31xl0 _29 (CKm); electronic polarizability, 

3. 8 A 3 (3.8xl0 -30 m 3 ). Total error, 5xl0' 3 and lxl0“ 2 ; 3 To 5 random number 
sets per point. 
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figure 2(a) is characteristic of the low velocity behavior observed in references 1 and 3. 

The Ctj curve becomes smoother at v = 10 , and the structure disappears entirely 

5 

for high velocity. For v = 2x10 centimeters per second the Ct, values drop off to 

o 

zero rapidly at b values beyond 6 A (0. 6 nm). 

The a values are plotted against initial ion energy in figure 3 and compared with 

values of the Langevin o, , permanent dipole ct_. and a cross sections. The ion 

energy is equal to twice the relative translational energy for this case. The cr values 
0 fifi 1 /9 ^ 

have a slope = ej , that is, somewhat steeper than the calculated e' 17 ^ for the 



Figure 3. - Comparison of numerically calculated capture cross section with 
various classical cross sections for methyl cyanide parent-ion collision. 


HCl + HC}. + system in reference 3. The symmetric top cr values are, however, much 
larger than the Langevin (simple polarizability) theory predicts. 

The integrated values of a (over all possible relative energies), designated as Q , 
are obtained from equation (2). 

The corresponding capture cross section assumed by the experimentalists is the 
integrated form of equation (4): 
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800 , 



Figure 4. - Comparison of integrated numerical capture 
cross section Q c with experimentally assumed capture 
cross section Q max and observed reaction cross 
section Q R for methyl cyanide parent-ion collision. 


The integrated numerical cross section, obtained by assuming that a c varies as the n^ 
power of the ion energy, is 


Q c 



"c< e l) d6 l = — 
e m 


l-n _ 1-n 

m i 

(1 -n) 


( 6 ) 


where a is a constant obtained from the assumed o e“ n dependence of the mono- 
energetic cross section. It should be noted that, for m < m^ (in all cases), the experi- 
mental values of ion energy will correspond to smaller values of relative energy. The 
numerical predictions have a slope nearly The experimentally observed reaction 

cross section is plotted for comparison with Q„ and <D „„ in figure 4. It should be 

noted that the term containing ef ' has been retained in computing Q„ and <3 „„ 

i c max 

but was not retained in the form of derived in reference 5. An effective reaction 
efficiency P~. of 0. 72 was assumed with the form of o in reference 5 to obtain 
Qj^; this precludes any simple interpretation of the disparity in results. Figure 5 shows 
cr c as a function of ion energy for the CHgCN + CgD^ collision. The inverse 
dependence on ion velocity for fixed reduced mass is identical to the results for the 
CH 3 CN parent ion collision shown in figure 3. Comparison of figures 3 and 5 indicates 
that a depends only on ion energy. From this fact one can infer that a depends 
roughly inversely on the mass for constant velocity v. 
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Case II: C^Dj2 + CH^CN 

As a result of the two dependences on v and m, the a values appear nearly 
independent of ion for the range of ion energies studied with CHgCN targets. This 
system was studied at the same ion velocities as case I, so the numerical points are 
merely shifted to higher because m^ is larger. The experimental reaction cross 
section Q R is plotted for comparison in figure 6. An effective reaction efficiency of 
0.67 was assumed in reference 5. Unlike the parent ion reaction, the values of Q R are 
greater than Q„ and Q„_„ at all energies. 

C lllciX 

Case III: Pseudo CH 3 CN + + CH3CN Collision 

The value of Ig/Ij chosen was 0. 24; all other constants were the same as for the 
case I. There were no significant differences in the C R curves between cases I and III. 
This indicates that the interaction is insensitive to the value of moment of inertia about 



Ion energy, eV 

Figure 7. - Comparison of monoenergetic capture cross 
section results for linear methyl cyanide (CH^CN + 

C H 3 C N + ) collision with results for hydrochloric acid. 
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the symmetry axis for symmetric tops. The corresponding cross section values ct c 
and Q do not differ from the case I results. 


Case IV: CH 3 CN + + "Linear CH3CN" 

This hypothetical linear molecule target was studied to investigate the parent ion 

collision in the limit where the rotational degree of freedom about the symmetry axis 

2 

vanishes. The against b curves differ only slightly from those of cases I and III 
at low velocity. However, there is a difference (of several percent) between cr c values 
in this region. The results are shown in figure 7 and are compared with similar results 
for the linear HC1 (ref. 3). The HC1 results have a smaller dipole moment and 
polarizability than CH 3 CN; a = 2. 6 A 3 (2. 6xl0 -30 m 3 ), \x = 1. 08 D. U. (3. 60xl0" 30 (C)(m)). 
The results indicate that the cross section does have the dipole moment dependence 
that the approximate theory predicts (refs. 1 to 3). The theory predicts values of o c 
nearly proportional to \±. Reference 3 also includes results for a hypothetical molecule 
(/ u and a equal to HC1 values) with an Ij value 100 times that of HC1. The results for 
linear CHgCN (Ij value ~ 40 times HC1) are in agreement with these predictions in that 
cr is nearly independent of the moment of inertia. 

V 


CONCLUDING REMARKS 

The interpretation of experimental results is complicated by recent experiments 
reported in reference 7 which indicate that there are competing reaction effects that have 
been ignored in obtaining ion molecule reaction cross sections. The numerical calcula- 
tions for symmetric top systems indicate only a slightly stronger dependence of capture 
cross section on ion energy than calculated for linear polar molecules. It is hoped that 
these numerical results might serve as a basis for interpreting future experiments on 
such reactions. However, more reliable absolute values of experimental reaction cross 
sections are necessary before comparison of theory and experiment can be confidently 
interpreted in terms of reaction efficiencies. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 21, 1968, 

120-26-03-04-22. 
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APPENDIX A 


EQUATIONS OF MOTION AND VALUES OF MOLECULAR CONSTANTS 


The Lagrangian energy expression of reference 1 for the ion-symmetric top molec- 
ular system which is depicted in figure 1 is 




+ 77 cos & 


+ — {-X sin £ sin 77 + Y sin £ cos 77 + Z cos £} + ii- 
r 3 2r 4 


(Al) 


where m, I, \i, a, and e are the reduced mass, moment of inertia, dipole moment, 
electronic polarizability, and electronic charge, respectively. 

The equations of motion of the system are 


d_ 

dt 




(A2) 


where q = X, Y, Z, 77, and )//. 

For the translational coordinates, equation (A2) reduces to the following equations: 


X = -Hi- {-r 2 sin | sin 77 - 3XF } - 2^X 
mr 3 mr® 

Y = -Hi {r 2 sin £ cos 77 - 3YF V } 

mr 5 mr 6 

Z = -Hi. {r 2 cos | - 3ZF } " 2ae2z 
5 v 6 


(A3) 


(A4) 


(A5) 


where the potential function 

F y = {-X sin £ sin 77 + Y sin £ cos 77 + Z cos £} 
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The corresponding rotational coordinate equations are 


i = i) 2 sin £ cos g {z sin 4 - X cos | sin ■q + Y cos 4 cos tj} 

v 3 


Io . . .2 

- — (\pn sin ij + t) sin 4 cos £) 


(A6) 


V 


)j-Me S j n 77 + X cos rf) - 21^1 cos | + I 2 [|(i// + 77 cos £)]] 

r^ 

Ij sin | 


(A 7 ) 


^ = 77^ sin (j - 77 cos £ 


(A8) 


It should be noted that the ion-dipole potential term in equation (Al) has a different 
form than the term used for the linear polar molecule in reference 3 . The dipole moment 
is assumed to lie along the symmetry axis of the top but the angle tj is equal to the con- 
ventional azimuthal angle plus 7r/2 (see fig. 2 ). This angle 77 locates the intersection 
of the plane normal to the symmetry axis in the x-y plane and not the projection of the 
symmetry axis in that plane. The form of the interaction potential term changes under 
the coordinate transformation (to the second Euler system), which was described in 
Details of the Calculation. 
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APPENDIX B 


VARIABLE STEP-SIZE ROUTINE 


There are two general classes of formulas for numerical integration of ordinary 
differential equations: the multistep (e. g. , predictor-corrector) and the single step or 
Runge-Kutta. Both types are discussed in reference 8. Multistep formulas are appealing 
because they yield a higher order of accuracy for the same number of evaluations of the 
differential equations. However, the Runge-Kutta formulas are more flexible in that 
they permit step-by-step adjustment of the integration interval. This adjustment 
feature is especially useful in a system where a reliable, independent indicator of inte- 
gration error is available. For a conservative interaction the deviation in total energy 
from its initial value serves as such a check because it can be easily calculated at the 
conclusion of each integration step. For this reason, Runge-Kutta formulas were used 
throughout the calculation. The error control and integration interval selection scheme 
used is essentially empirical; it works well, although it is without a rigorous mathemati- 
cal basis. The error is considered through Rj., the ratio of the current total energy (at 
time t) to the initial energy of the system. The error E t to be tolerated in any step 

is simply the absolute value of the change in R. during the step. A value E™, of the 
-4 9 t A 

total error (10 < E^, < 10“ ) to be tolerated over the total estimated trajectory time 

T is selected. It is then assumed that the error accumulates monotonically, that is, 

no cancellation of error. The acceptable error for one integration step of length h at 

time t (with current energy ratio Rj.) is 


E t =h 


(E T - JRj. - 1.00 . . . |) 
__ 


(Bl) 


If the energy ratio at the end of this step is R^, then the actual error accrued during 
this step is 


E a =K-Rtl < B2 > 

If E a > CjE t , where 1. 0 < Cj < 4. 0, then the step is not accepted. 

After both successful and unsuccessful steps, one must determine the size h’ of 
the next step. If the step is unsuccessful, the step should be repeated with h' < h; and 
if the step is successful, the next step should be performed with h' > h. For a Runge- 
Kutta formula of order m applied to the single differential equation, 

y’ = f(t, y) (B3) 
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the truncation error is given by 


E m = Ah m+1 + 0(h m+2 ) (B4) 

where A is a function of the formula used and of f(t,y). 

Assuming that E a = Ah m+ * (not theoretically justified), then the value of A can be 
estimated. With A known, the assumption 


e; = Ah' m+1 =. Ea h ,m+1 (B5) 

h m+l 

where represents the error to be accepted in the next step, yields 

E , l/(m + 1) 

h' = — • h (B6) 

E a 

A limitation on h* which has proved useful and, in fact, necessary for E = 0, is 

a 

h' < c 2 h (B7) 

20 

where 2. 0 < c 2 < 10. 0. In addition, any case is terminated if h* < 10 seconds. 

This situation is rarely encountered if realistic values of E T and T have been assigned. 

This procedure for selecting integration intervals satisfies the previously stated 
qualitative criteria. If an unrealistic estimate is given for the initial step, the step size 
settles down to a sensible value after two or three integrations (or failures, if the esti- 
mate is too large). The scheme has yielded more accurate results and significantly 
shorter computer run times than fixed-step size routines. Calculations were carried 
out using Runge-Kutta formulas of order 2 (requiring 2 evaluations of the differential 
equations per step) and 3 (3 evaluations) both described in reference 9. Calculations 
were mainly done with orders 4 (4 evaluations - the classic Runge-Kutta method) and 5 
(6 evaluations) which are described in references 8 and 10, respectively. Formulas for 
orders 6 (7 evaluations), 7 (9 evaluations) and 8(12 evaluations) are described in refer- 
ences 11 and 12. 
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